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摘 要 :采用 精确 解析 求解 方法 研究 石墨 烯 片 (GPL) 增 强 功 能 梯度 梁 的 非 线 性 屈曲 及 后 屈曲 力学 
行为 。 研 究 中 假定 GPL 均匀 地 分 布 在 梁 的 每 一 层 , 但 其 重量 分 数 沿 厚度 方向 呈 梯 度 或 均匀 变化 , 导 
致 结构 的 物性 参数 沿 厚 度 方 向 连续 梯度 变化 ,并 根据 Halpin-Tsai 力学 模型 和 混合 率 法 则 得 到 材料 
,的 等 效 物性 参数 。 同 时 基于 Euler-Bernoulli 梁 理 论 和 Von-Karman 非 线 性 应 变 位 移 关 系 由 能 量变 分 
法 获得 石墨 烯 片 增强 功能 梯度 梁 的 非 线 性 控制 方程 ,对 此 方程 进行 精确 解析 求解 获得 临界 载荷 、 慑 
曲 模 态 以 及 后 属 曲 变形 的 解析 表达 式 。 并 进行 参数 研究 分 析 了 GPL 的 重量 分 数 、 几 何 形状 及 尺寸 
和 分 布 方 式 对 该 复合 材料 梁 非 线 性 屈曲 和 后 屈曲 的 影响 。 
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Exact solutions for buckling and post-buckling of functionally 
grade graphene reinforced composite beams 
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2. Liuhuan Transmission ( Xian) Technology Co.,Ltd,710000 Xi’an, China) 


Abstract : In this paper, nonlinear buckling and post-buckling behaviors of functionally graded beams rein- 
forced by graphene platelets (GPL) are studied by exact analytical method. It is assumed that GPL is uni- 
formly distributed in each layer of the beam, and its weight fraction varies gradient along the thickness di- 
rection , which leads to continuous and gradient changes of the materials properties along the thickness. Ac- 
cording to Halpin-Tsai mechanical model and mixing rate rule ,the equivalent materials properties are ob- 
tained. Simultaneously , based on Euler Bernoulli beam theory and Von-Kérman nonlinear strain-displace- 
ment relationship, nonlinear governing equations of the graphene platelets reinforced functionally graded 
beam are obtained by energy variational method. And the analytical expressions of critical buckling loads, 
buckling modes and post-buckling deformations are obtained through exactly and analytically solve the 


nonlinear governing differential equations. The effects of the weight fraction, geometrical parameters and 
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distributions of GPL on the nonlinear buckling and post-buckling of the composite beam are analyzed by 


parameters research. 
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石墨 烯 独 特 的 二 维 结构 赋予 它 优异 的 力学 性 
能 ,使 它 成 为 非常 理想 的 增强 材料 。 将 少量 的 石墨 
烯 或 其 衍生 物 作为 增强 体 加 入 基体 成 为 石墨 烯 增强 
复合 材料 中 , 比 传 统 碳纤维 和 碳 纳米 管 更 能 显著 提 
高 结构 的 强度 和 刚度 。 将 功能 梯度 材料 卓越 
的 设计 理念 推广 到 石 暴 烯 增强 复合 材料 中 ,也 就 是 
将 石墨 烯 片 (graphene platelets, GPL) BAK AM /) ERE 
为 增强 体 方向 随机 且 均 匀 地 散布 在 每 层 基 体 中 , 石 
墨 烯 含量 在 基体 厚度 方向 上 逐 层 呈 梯 度 或 均匀 分 
布 ,形成 新 一 代 的 先进 复合 材料 , 即 功能 梯度 石墨 炳 
增强 复合 材 Æ} U (functionally graded graphene plate- 
let& reinforced composite, FG-GPLRC ) , FG-GPLRC 
也 吉 有 非凡 的 力学 性 能 ,主要 取决 于 石墨 烯 ,可 以 通 
过 会 理 设计 石墨 烯 片 的 分 布 方式 、 重 量 分 数 和 几何 
参数 得 到 最 优 的 整体 力学 性 能 。 
QOFC-GPLRC 结构 的 力学 性 能 研究 近年 来 已 成 为 
国体 力学 中 一 个 新 的 研究 领域 。 在 线性 稳定 性 方 
HL KITIPORNCHAI a RH EZX GPL 增强 功 
能 梯度 多 孔 梁 的 自由 振动 和 弹性 届 曲 进行 了 研究 。 
VORES djab FAIRS H N ZA REENE 
层 在 墨 烯 增强 功能 梯度 梁 的 动力 稳定 性 进行 了 研 
FEE YANG 等 'q 基 于 能 量 法 研究 了 石墨 烯 增强 纳米 
复 叭 材料 拱 在 点 载荷 作用 下 的 动态 屈曲 。LU 等 中 
基于 一 阶 剪 切 变形 板 理论 和 最 小 势能 原理 研究 了 石 
墨 医 增强 功能 梯度 压 电 板 的 届 曲 特性 。KIANIW 对 
热 环 境 中 受 外 部 压力 作用 下 功能 梯度 石墨 烯 增强 锥 
过 的 屈曲 进行 了 研究 。TAM 等 ' 对 含有 开放 边缘 
裂纹 的 石墨 烯 增强 功能 梯度 复合 材料 梁 的 届 曲 进行 
了 研究 。SHAHGCHOLIAN 等 (中 基于 一 阶 前 切 变形 
理论 对 石墨 烯 增 强 多 孔 纳米 复合 材料 圆柱 壳 的 屈曲 
进行 了 研究 。PARASHAR 等 "| 采用 有 限 元 法 对 受 
轴 问 载荷 作用 下 的 石墨 烯 增强 环 氧 复 合板 的 届 曲 进 
行 了 研究 。 通 过 研究 表明 当 石 墨 烯 片 的 体积 分 数 是 
6% 时 ,该 结构 的 届 曲 强度 提高 26% 。EBRAHIMI 
等 (中 研究 了 热 环境 中 GPL 增强 聚合 物 复合 材料 板 
在 黏 弹性 地 基 上 的 屈曲 特性 。 同 时 考虑 了 GPL 的 4 
种 分 布 模式 ,通过 Hamilton 原理 建立 了 GPL 增强 板 
的 控制 方程 ,用 Navier 解法 对 该 问题 进行 了 分 析 求 
解 。SUN 4°19) 基于 Donnell 壳 理 论 对 石墨 烯 增强 功 
能 梯度 复合 材料 圆柱 壳 的 扭转 屈曲 进行 了 研究 ,并 


由 此 得 到 了 FG-GPLRC 壳 体 的 临界 屈曲 载荷 和 相应 
的 屈曲 模 态 。 

以 上 研究 仅 限于 线性 问题 。 由 于 FG-GPLRC 结 
构 的 高 强度 .高 韧性 使 其 能 承受 较 大 变形 ,对 该 类 结 
构 进行 非 线性 大 变形 分 析 远 比 线性 小 变形 分 析 更 加 
有 意义 。 在 非 线性 稳定 性 的 研究 方面 ,YANG 等 中 
在 考虑 弹性 地 基 影 响 的 同时 通过 微分 求 积 法 对 石墨 
烯 增强 功能 梯度 复合 材料 梁 的 屈曲 及 后 屈曲 问题 进 
行 了 研究 。CHEN 4U php Fy ELBE RE ARIE OTE T 
沫 金属 基 GPL 增强 多 层 功 能 梯度 多 孔 纳米 复合 材 
料 梁 的 后 屈曲 行为 。SHEN 等 "研究 了 功能 梯度 石 
墨 烯 增强 复合 材料 层 合 梁 在 弹性 地 基 上 的 热 后 屈曲 
及 非 线 性 弯曲 。WANG 26") 对 GPL 增强 压 电 复合 
材料 梁 的 屈曲 和 后 屈曲 行为 进行 了 分 析 , 考虑 剪 切 
变形 和 非 线性 应 变 -位 移 关 系 ,利用 虚 功 原理 建立 了 
该 梁 的 控制 微分 方程 ,采用 微分 求 积 法 定量 地 研究 
影响 届 曲 性 能 的 因素 。ANSARI 等 :的 基于 一 阶 剪 
切 变形 过 理论 和 Von-Karman 非 线性 几何 关系 研究 
了 功能 梯度 石墨 烯 片 增强 聚合 物 纳米 复合 材料 圆柱 
壳 的 后 届 曲 ,并 利用 Fourier 级 数 和 微分 求 积 数 法 提 
出 了 半 解 析 解 的 方法 。SONG 等 i' 在 考虑 弹性 地 基 
和 边缘 裂纹 影响 的 同时 基于 一 阶 剪 切 变形 梁 理 论 ， 
研究 了 石墨 烯 复合 材料 多 层 梁 的 热 届 曲 和 后 屈曲 。 
SAHMANI 等 ”1 基于 非 局 部 应 变 梯度 理论 研究 了 轴 
向 载荷 作用 下 GPL 增强 多 层 功 能 梯度 复合 纳米 板 
的 非 线性 屈曲 及 后 屈曲 。 

目前 石墨 烯 增强 功能 梯度 材料 结构 非 线 性 力学 
行为 的 已 有 研究 还 比较 少 ,尤其 对 该 结构 非 线性 稳 
定性 问题 的 非 线性 控制 微分 方程 进行 精确 解析 求解 
方面 鲜 有 文献 报道 。 本 研究 在 考虑 非 线性 大 变形 的 
情况 下 ,基于 Halpin-Tsai 力学 模型 2 ,用 解析 法 精 
确 求解 石墨 烯 增强 功能 梯度 粱 的 非 线性 稳定 性 控制 
方程 ,获得 临界 屈曲 载荷 . 届 曲 模 态 和 后 屈曲 变形 的 
表达 式 ,并 分 析 GPL 的 重量 分 数 .不 同 的 分 布 模式 、 
几何 形状 和 尺寸 等 对 临界 屈曲 载荷 的 影响 。 


1 数学 模型 


1.1 材料 物性 参数 
考虑 图 1 ERKA LIEN b AN h 的 石墨 烯 增 
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成 ,并 且 每 层 中 的 石墨 烯 为 均匀 分 布 ,但 不 同 层 中 石 
墨 炳 含量 不 同 。 本 研究 考虑 各 层 中 石墨 烦 的 含量 具 
有 U-GPLRC ,X-GPLRC 和 0-GPLRC 3 种 不 同 分 布 
模式 如 图 1 所 示 , 图 中 颜色 越 深 则 石墨 烯 浓度 越 大 。 
U-GPLRC 表示 GPL 在 每 一 层 均 匀 分 布 且 浓 度 相同 ; 
X-GPLRC 表示 GPL 在 上 层 和 下 层 具有 较 高 的 浓度 
而 在 中 间 浓 度 最 低 ;0-GPLRC 表示 GPL 在 上 、 下 层 
有 具 有 较 低 的 浓度 而 中 间 层 具有 较 高 的 浓度 。 这 3 种 
分 布 模式 的 第 下层 GPL 体积 分 数 Vem 表达 分 别 为 
U-GPLRC: Vin, = Veo. » 

X-GPLRC; Vi, = 2V [2k —N, -1|/NL. 
O-GPLRC; Von, = Vén - |2k - N, — 197M. 


IF: k = 1,2,3…N ，N 是 梁 的 总 层 数 ; Voy, 是 
总 的 体积 分 数 , 可 表示 为 
© p Sori 
O Vepr = 
© gcpL + Cpcp./pm? 1 = gcpL) 
HD: 8 cpL 是 石墨 烯 的 质量 分 数 ; Pu All Pop. 分 别 为 
人 墨 烯 的 密度 。 


U-GPLRC X-GPLRC 


O-GPLRC 


图 1 3 种 不 同 的 GPL 分 布 模式 
Fig.1 Three different GPL distribution patterns 


FG-GPLRC 结构 第 层 的 弹性 模 量 £5” 基于 目 
前 广泛 采用 的 Halpin-Tsai 微观 力学 模型 给 出 "3 , 即 


EP = ( 31 + ENVE 51 Ua 
c cin + Ck M 
8 1- NLV cr 8 1- Mw Verr 


chinaXiv:: 


(1) 
> (Eom/Ew -1 _ CE op/Ey? ~ 1 
oe m= CEcp /Ey) + éL’ AS CE op /Ey) + Ew : 
Ey 和 E cr 分 别 为 基体 和 石墨 烯 的 弹性 模 量 ， ÉL All 
Ey 表示 GPL 的 几何 形状 尺寸 参数 ,具体 为 
= lop. _ 2/ Scr. 
é = (2), by 2( a) 
hop 和 bop, 分 别 表 示 GPL 的 长 度 宽度 和 


hep. 

式 中 /cr 
厚度 。 

FG-GPLRC 的 泊 松 比 可 由 混合 定律 给 出 


J 
a fE m hua mR NAO ERT 591 


Ck) ck k 
Vo = Vor Ver, + UML — Vege) (2) 


式 中 Vy 和 VepL 分 别 为 基体 材料 和 石墨 烯 的 泊 松 比 。 
12 梁 的 控制 方程 


根据 Euler-Bernoulli 粱 理论 , 同时 考虑 Von- 
Kármán 非 线性 应 变 位 移 关系 ,可 以 得 到 梁 的 轴 向 应 
变 表达 式 为 
e, = +t (%) + 2K (3) 
式 中 ; x 和 = 分 别 为 轴 向 和 横向 坐标 ; A w 分 别 为 
铀 向 位 移 和 找 度 ; x = - CF 为 梁 中 性 层 的 曲率 。 


考虑 线 弹 性 范围 内 的 变形 ,由 弹性 本 构 关 系 可 
1345 k JÆ FG-GPLRC WIEN A! 
Ch _ Ec » du 1 dw 2 dw 
| 
将 其 沿 厚 度 方 向 进行 积分 ,可 获得 粱 横 截 面 上 
的 轴 力 N, MEE M, 表达 式 , 即 


N, = foda = 4, [42+ (4) ] Ga) 
M sf dA =- D, dw 
= Fae =-= DV die (Sb) 


式 中 : 4 为 梁 横 截面 面积 ; 4， 和 Di 分 别 表示 梁 的 
拉 压 刚度 和 老 曲 刚度 ,通过 积分 计算 可 得 


N, ck 
Ec h po 
{A,, Dy} of mi :| 
iiaii pa Lao? N? 3 


采用 最 小 势能 原理 67 = 8U + V) = 0 可 以 推 
导出 梁 的 控制 方程 , 式 中 开 、V 和 了 分 别 表示 系统 
的 总 势能 ,应变 能 和 外 力 势 能 , 8 表示 变 分 。 

FG-GPLRC 梁 的 应 变 能 及 其 一 阶 变 分 为 


1 1y i 
U = ratat = 和 Ik N 8e? + M.8x)dx 


spet = Sg F(t)” 。 考 虑 梁 上 作用 给 向 压缩 


V = - [4p( 22) ax see P Pes( te) ax 
小 势能 原理 得 到 梁 的 控制 方程 为 


dN, 
ae = 0 (6) 


ae 
(et 
地 


dw Ëw _ 
-Dı (fee 0 (7) 
X} FG-GPLRC 深 的 屈曲 及 过 屈曲 问题 的 研究 ， 
即 对 方程 组 (6) 和 (7) 进 行 求解 ,还 需 同 时 考虑 边界 
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条 件 。 在 本 研究 中 考虑 梁 的 两 端 均 是 夹 紧 边 界 , 相 
应 的 边界 条 件 为 


00 (8) 


1.3 精确 解析 求解 


对 给 定 边 界 条 件 的 非 线性 控制 微分 方程 边 值 问 
题 进行 精确 解析 求解 。 首 先 将 式 (6)、 式 (7) 采 用 位 
移 形式 表示 为 


reele 9) 
Bu Ge + (aul Ge ty (Se) 1 Pfs = 
(10) 
“PGE (9) X x 进行 一 次 积分 可 得 
O A, [4 + 3 (%) | = B, (11) 
加 二 次 积分 为 
= Ay [u + [> (4) a] 2RaaBs . (ia) 


Dou T oR AR B, ALB, ,将 边界 条 件 (8) 代入 式 
(12) 可 得 


.二 再 将 B, 所 得 值 代入 式 (11) 得 到 
~~ du 1 dw 2 E An 1 dw 2 
Solat rla) l= aha) & 
ER ERRAR (10 ) JR B a Ae 8 BE 
RHIA, BI 
dfw p Ay : dw : dw = 
Die aria) “ae TA 
为 了 使 所 求 的 解 具有 一 般 性 ,引入 无 量 纲 量 , 即 
w Anh? PP D, 
er Dy “Da Dg 


通过 代入 无 量 纲 量 , 可 获得 石墨 烯 增 强 功 能 梯 
度 梁 拨 度 的 无 量 纲 控制 方程 , 即 
dw F rydW vldW_ 
| | =0 (14) 


dX dX’ 
夹 紧 端的 无 量 纲 形式 边界 条 件 为 
dW 
X=0,1h,W = 7 =0 


为 了 使 所 求 的 方程 形式 上 更 加 简洁 , 令 


parE fafa a 
则 式 (14) 转 化 为 
dd -0 (16) 
方程 (16) 的 通 解 为 
W(X) = Cicos(AX)+ Cysin(,X) + CX + C, 
(17) 


AP C, 0, CG 和 C, 表示 积分 常数 。 进 一 步 将 上 
述 边界 条 件 代 入 式 (17) 中 有 
Wo) = C, + C, = 0 
W(1) = C,cosA + CysinA + C, +C, = 0 
(18b) 
(18c) 


(18a) 


Wo = CA+C =0 
W (1) =- CAsinA + C,AcosA + C; = 0 


(18d) 

令 C, = C , 联 立 求解 式 (18a) ~ (18d) 可 得 

B cosà - 1 _ (1 - cosà)À 
r= a Os ee sind 一 从 fks 


此 时 ,精确 解析 求解 获得 了 在 该 边界 条 件 下 石 
墨 烯 增强 功能 梯度 梁 届 曲 问 题解 的 形式 , 即 
cosA - 1 


sind 一 从 


WX) = cl - cos( AX) + sin (AX + 


ACE = eh) 
WP C 为 与 轴 辐 载荷 已 有 关 的 常数 。 再 将 式 (19 ) 
代入 式 (15 ) 可 得 
pan ee (20) 
JF NA 
替换 式 (18d) 中 的 Ci, Ca AC, 可 得 边界 条 件 
为 两 端 夹 紧 时 ,石墨 烯 增强 功能 梯度 梁 届 曲 的 特征 
方程 , 即 


2 -2coshA -AsinA = 0 (21) 

由 特征 方程 (21) 可 以 求 得 边界 条 件 为 两 端 夹 
紧 时 的 无 量 纲 特 征 值 为 

A, = 6.283 , A, = 8.987 ,A, = 12.566, A, = 
15.451 , A; = 18.850 ,--- 

采用 归 一 化 方法 令 C = 1 ,将 特征 值 A, 代入 式 
(19) 可 得 石墨 烯 增强 功能 梯度 梁 的 第 n 阶 届 曲 
模 态 。 
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(cosA, — 1)(sinA,X - AX) 
sind, — A 


W, = (1 - cosd,X) + 


(22) 

由 最 小 特征 值 4，= 6.283 , 即 为 27 ,将 其 代入 
到 式 (15) ,可 得 两 端 夹 紧 时 FG-CPLRC 梁 的 届 曲 临 
界 载荷 的 解析 表达 式 , 即 


D 
P =4,/ = 
a T D, 


m2 è Em 2,3 z 
RPD, = | rz dz HAIER RIE, 


-12 1 -— vy 

对 齐 次 方程 (16) 进行 求解 可 得 两 端 夹 紧 
FC-GPLRC 深 的 临界 屈曲 载荷 对 应 的 后 屈曲 挠 度 表 
达 式 为 


(x) 2 [fo Po 1[ ] 
) = JF dg D, 1 - cos(2TX) 


(23) 


树脂 作为 基体 材料 ,其 弹性 模 量 为 3 GPa 、 密 度 为 
12G em 、 泊 松 比 为 0.34) 。 综 合 考虑 制造 成 本 


1.00 TPa 密度 为 1.062 g/ cm? 、 泊 松 比 为 0.186。 如 
无 疆 别 说 明 ,石墨 烯 的 质量 分 数 为 0.3% ,几何 尺寸 为 
E 25 pm bon, = 1.5 pm Al hen, =1.5hm >。 

为 了 证 明 本 研究 的 理论 推导 及 求解 过 程 的 正确 
性 ,将 本 研究 的 结果 与 文献 114 ] 中 使 用 微分 求 积 法 
所 得 的 相应 结果 作 比 较 。 引 入 文献 [14] 中 的 无 量 


纲 载荷 形式 : P = 1 ,其 


表示 纯 基 体 梁 的 刚度 系数 。 计 算得 到 两 端 夹 紧 边 界 
条 件 下 , 当 GPL 的 重量 分 数 分 别 为 0.1% ALO. 3% 时 
梁 的 无 量 纲 届 曲 临界 载荷 的 数值 结果 ,并 将 该 临界 
载荷 与 文献 14] 中 使 用 微分 求 积 法 所 获得 的 相应 
结果 同时 列 入 表 1 中 。 由 表 1 可见, 用 本 人 研究 中 的 
解析 求解 方法 所 得 的 结果 与 文献 [ 14 ] 中 数值 方法 
所 求 的 结果 极其 接近 , 故 证 明了 本 研究 理论 推导 及 
解析 求解 过 程 的 正确 性 。 
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表 1 石墨 烯 增强 功能 梯度 梁 的 无 量 纲 屈曲 临界 载荷 比较 


Tab.1 Comparison of dimensionless critical buckling 


loads for the FG-GPLRC beams 


Zer./%  U-GPLRC X-GPLRC 0-GPLRC 

本 研究 0.1 0.038 4 0.041 0 0.033 2 
文献 [14] 0.1 0.0391 0.0433 0.0349 
本 研究 0.3 0.0577 0.0703 0.044 1 
文献 [14] 0.3 0.058 6 0.0709 0.045 8 


FACILE A shad ah IE BES RIAIT GT De HH 
模 态 (图 2) ,由 图 2 可 见 该 复合 材料 梁 的 届 曲 构 型 与 均 
匀 材 料 或 者 一 般 功能 梯度 材料 梁 的 届 曲 构 型 相同 。 


1.0 Bee 
Fas 


en. 
POENA AS 
Eg \ 
4% 


rf 2 
7 


\ 


0.0 0.2 0.4 0.6 0.8 1.0 


图 2 石墨 烯 增强 功能 梯度 梁 的 前 4 阶 届 曲 模 态 
Fig.2 The first four order buckling modes 
of the FG-GPLRC beam 
其 次 分 析 FG-GPLRC 梁 的 临界 届 曲 载荷 及 其 影 
响 因素 。 图 3 给 出 了 不 同 石墨 烯 重量 分 数 梁 的 临界 
屈曲 载荷 ,其 中 横 坐 标 和 纵 坐 标 分 别 表 示 GPL 的 宽 
FEE EAE KS Il AJE HH ER TT o 


400 


X-GPLRC 
350b Lon /Don=1, Uh=10 


Ban!hen. 
图 3 A2 he Ce ay BOOT lh FL Hh ee Hi Be n 


Fig.3 Effect of graphene mass fraction 


on the critical buckling loads 
由 图 3 可 见 当 GPL 在 基体 中 的 占 比 越 高 时 , 梁 
的 临界 屈曲 载荷 也 越 大 。 尤 其 当 石 墨 烯 片 的 宽厚 比 
较 大 时 ,增强 效果 更 加 明显 ,石墨 烯 片 宽厚 比 为 
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10° .含量 仅 为 0.3% 时 的 FG-GPLRC 梁 的 临界 届 曲 
载 向 是 纯 基 体 梁 的 4.2 倍 。 

为 了 更 进一步 考察 石墨 烯 片 的 几何 尺寸 对 梁 届 
曲 的 影响 ,图 4 画 出 了 梁 的 临界 届 曲 载荷 随 GPL 的 
几何 参数 变化 的 关系 曲线 ,其 中 图 4(a) 和 图 4 Cb) 
的 横 坐 标 分 别 表示 GPL 的 宽厚 比 和 长 宽 比 , 纵 坐 标 
为 梁 的 临界 届 曲 载荷 。 


400 


X-GPLRC 
Sam=0.3%, Up=10 


350 


-一 一 = 


om) 
ot 


250 
200 
b 
> 
150 
A o ™ P i 
© Dis. Mess: 
= (a) 宽厚 比 的 影响 
© 280 
© BonW03%, bean/hon100  - ------------ 
© 2% 
N ee 
= 0 dt on gee 
N 
E u 
a ---- X-GPLRC 
; = 一 = O-GPLRC 
>< 
200% 
二 0 5 T T a 
f = lonl bon 
PT (b) 长 宽 比 的 影响 


Qaa 石墨 烯 片 的 几何 尺寸 对 临界 届 昌 载荷 的 影响 
Fig.4 The influence of the geometry size of 
CPL on the critical buckling loads 

在 石墨 烯 总 含量 和 石墨 烯 的 宽度 不 变 的 情况 
下 , 较 大 的 长 宽 比 意味 着 石墨 烯 片 具有 和 较 大 的 表面 
积 , 较 大 的 宽厚 比 意味 着 在 石墨 燃 中 所 含 的 单 层 石 
墨 烯 片 较 少 '% 。 因 此 ,从 图 4(a) 中 可 以 发 现 随 着 宽 
厚 比 的 增 大 ,临界 屈曲 载荷 也 随 之 增 大 ,说明 含有 较 
少 的 单 层 石墨 烯 片 时 增强 效果 更 加 明显 。 但 当 宽 厚 
比 大 于 10 时 这 种 变化 不 太 明 显 , 即 宽厚 比 大 于 10° 
时 ,再 无 法 通过 增 大 GPL 的 宽厚 比 来 有 效 提高 梁 的 
强度 。 从 图 4(b) 中 看 到 随 着 长 宽 比 的 增 大 ,临界 届 
曲 载荷 也 随 之 增 大 ,说 明 具 有 较 大 表面 积 的 石墨 烯 
片 对 于 功能 梯度 梁 的 强度 具有 明显 的 增强 效果 。 但 
当 lop/bem > 4 后 对 临界 屈曲 载荷 的 影响 变 得 不 太 
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明显 , 即 lw/bcwm > 4 时 ,也 无 法 通过 增 大 GPL 的 长 
宽 比 来 提高 增强 效果 。 由 图 4(b) 也 发 现 ,不 同 的 
GPL 分 布 形式 对 临界 屈曲 载荷 的 增强 效果 不 同 ,X- 
GPLRC 模式 能 更 有 效 地 提高 梁 的 临界 载 答 ,增强 效 
果 最 好 , 即 在 梁 的 上 下 表面 分 散 由 较 少 单 层 石 墨 烯 
片 组 成 的 石墨 烯 是 提高 粱 的 临界 载荷 的 最 有 效 
方法 。 

最 后 分 析 FG-GPLRC 梁 的 后 届 曲 响应 及 其 影响 
因素 。 图 5 和 图 6 分 别 绘 出 了 GPL 的 质量 分 数 和 分 
布 形式 不 同时 的 石墨 烯 增强 功能 梯度 梁 的 后 届 曲 路 
径 , 其 中 横 坐 标 和 纵 坐 标 分 别 表示 梁 的 轴 向 载荷 和 
后 届 曲 最 大 接 度 。 由 图 5 可 见 , 随 着 GPL 质量 分 数 
的 增 大 , 梁 的 后 届 曲 承载 能 力也 随 之 显著 增 大 。 由 
图 6 可 见 ,3 种 分 布 模式 下 的 FG-GPLRC 梁 都 具有 
稳定 的 后 届 曲 路 径 , 并 且 X-GPLRC 分 布 模式 下 梁 具 
有 更 高 的 后 届 曲 承载 能 力 , 其 次 是 U-GPLRC 分 布 ， 
最 后 是 O-GPLRC 分 布 。 也 就 是 说 X-GPLRC 分 布 模 
式 在 上 表面 和 下 表面 具有 更 高 的 GPL 含量 , 梁 的 弯 
曲 刚 度 最 高 ,所 以 梁 具有 更 强 的 承载 能 
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Fig.5 Effect of the graphene mass fraction 
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K6 3 种 GPL 分 布 模式 梁 的 后 届 曲 路 径 
Fig.6 Post-buckling paths of the beams 
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图 7 和 图 8 表示 了 GPL 的 几何 尺寸 对 后 屈曲 路 
径 的 影响 ,图 中 不 同 曲线 分 别 代 表 GPL 的 长 宽 比 和 
宽厚 比 不 同 梁 的 后 屈曲 路 径 。 由 图 可 见 , 随 着 GPL 


张 靖 华 ， 等 ; 功能 梯度 石墨 烯 增强 


nay iAd HE: 
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由 其 变化 曲线 可 知 不 同 GPL 4AE PR IN AY) Je JH 
曲 变形 有 一 定 的 差异 ,这 是 因为 不 同 分 布 模式 梁 的 
强度 不 同 , 相 比 其 他 分 布 模式 ,X-GPLRC 分 布 模式 


的 长 宽 比 和 宽厚 比 的 增 大 , 深 的 后 屈曲 最 大 挠 度 逐 
渐 减 小 ,具有 更 高 的 后 屈曲 承载 能 
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Fig.7 Effect of geometry size of graphene platelets 
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on the post-buckling paths 
图 9 绘 出 了 不 同 分 布 形 式 的 石墨 烯 增强 功能 梯 
度 梁 的 后 届 曲 构 型 图 ,其 中 给 定 无 量 纲 载 寄 P=250。 
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08 IEE O-GPLRC re s besi /ham=100 


图 9 FG-GPLRC 梁 的 后 届 曲 构 型 
Fig.9 Post-buckling configurations of the FG-GPLRC beams 
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本 研究 通过 精确 解析 求解 非 线 性 微分 方程 研究 
了 石墨 烯 增强 功能 梯度 梁 的 届 曲 及 后 届 曲 特性 , 获 
得 了 屈 曲 临界 载荷 `. 届 曲 模 态 及 后 屈曲 挠 度 的 解析 
表达 式 。 并 通过 解析 表达 式 的 数值 算 例 探讨 了 不 同 
的 GPL 分 布 模式 、 几 何 尺 寸 和 重量 分 数 等 参数 对 石 
墨 烯 增强 功能 梯度 梁 的 非 线 性 届 曲 及 后 届 曲 的 影 
响 , 获 得 如 下 主要 结论 。 

1) 加 入 极 少 量 的 GPL 可 以 显著 提高 功能 梯度 
梁 的 临界 屈曲 载荷 和 后 屈曲 承载 力 , 即 对 梁 的 强度 
和 刚度 具有 明显 的 增强 效应 。 

2) 当 GPL 具有 较 大 的 表面 积 和 较 少 的 单 层 石 
黑 烯 片 时 ,能 提高 对 梁 的 增强 效应 ,其 临界 载荷 增 
大 ,后 屈曲 挠 度 减 小 。 

3) 不 同 的 GPL 分 布 形 式 对 梁 的 增强 效果 不 同 ， 
当 石 黑 烯 以 X-GPLRC 分 布 时 , 即 梁 的 上 下 表面 具有 
较 多 的 GPL, 梁 具有 更 高 的 承载 能 力 以 及 临界 届 曲 
BAW o 
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